Isospin Dynamics in Heavy Ion Collisions 



O 
O 

(N 



> 
in 
m 
o 

o 

\o 
o 



o 

=5 



X 



M.Di Toro 
Laboratori Nazionali del Sud INFN, 
Via S.Sofia 62, 95123 Catania, Italy, 
and Physics-Astronomy Dept., University of Catania, 
E-mail: ditoro@lns.infn.it 

Some isospin dynamics results in heavy ion collisions from low to relativistic energies obtained 
through transport approaches, largely inspired by David M. Brink, are rewieved. At very low ener- 
gies, just above the Coulomb barrier, the stimulating implications of the prompt dipole radiation in 
dissipative collisions of ions with large isospin asymmetries are discussed. We pass then to the very 
rich phenomenology of isospin effects on heavy ion reactions at intermediate energies (few A GeV 
range). We show that it can allow a "direct" study of the covariant structure of the isovector in- 
teraction in the hadron medium. We work within a relativistic transport frame, beyond a cascade 
picture, consistently derived from effective Lagrangians, where isospin effects are accounted for in 
the mean field and collision terms. Rather sensitive observables are proposed from collective flows 
("differential" flows) and from pion/kaon production (7r~/7r + , K°/K + yields). For the latter point 
relevant non-equilibrium effects are stressed. The possibility of the transition to a mixed hadron- 
quark phase, at high baryon and isospin density, is finally suggested. Some signatures could come 
from an expected "neutron trapping" effect. The importance of violent collision experiments with 
radioactive beams, from few AMeV to few AGeV, is stressed. 



INTRODUCTION: THE DAVID LEGACY 



It is a great pleasure and honor to contribute to this 
special day devoted to celebrate the 75th anniversary of 
David Brink. After many years of acquaintance and col- 
laboration I would like to state a few points I mostly 
got from him: 1.) The Mean Field "moves" the world, 
particularly true in nuclear dynamics, as I will show in 
examples discussed here; 2.) The role of relativity in 
nuclear structure and reactions; 3.) Play always a deep 
attention to the suggestions of young people; 4.) Fol- 
low a "sensible" behaviour towards the political choices, 
i.e. try to change something only when you can count 
on something better. Here of course I will focus on the 
physics part showing a series of results broadly inspired 
by the David ideas. 

In the last years the isospin dynamics has gained a lot 
of interest, as we can see from the development of new 
heavy ion facilities (radioactive beams), for the unique 
possibilities of probing the isovector in medium interac- 
tion far from saturation, relevant for the structure of un- 
stable elements as_well as for nuclear astrophysics see the 
recent reviews 111 



Here I will show some selected results of the mean 
field transport approach in a wide energy range, from 
few AMeV to few AGeV, in non-relativistic and rela- 
tivistic frames. At low energies I will discuss the isospin 
equilibration in dissipative collisions, fusion and deep- 
inelastic, through a related observable, the Prompt Col- 
lective Dipole Radiation. At high energies I will shortly 
present isospin effects on collective flows, on particle pro- 
duction and finally on the transition to a mixed hadron- 
quark phase at high baryon density. 



II. THE PROMPT DIPOLE 7-RAY EMISSION 

The possibility of an entrance channel bremsstrahlung 
dipole radiation due to an initial different N/Z distri- 
bution was suggested at the beginning of the nineties 
0, 3, largely inspired by David discussions. At that 
time a large debate was present on the disappearing of 
Hot Giant Dipole Resonances in fusion reactions. David 
was suggesting the simple argument that a GDR needs 
time to be built in a hot compound nucleus, meanwhile 
the system will cool down by neutron emission and the 
GDR photons will show up at lower temperature. The 
natural consequence suggested in |3J was that we would 
expect a new dipole emission, in addition to the statisti- 
cal one, if some pre-compound collective dipole mode is 
present. After several experimental evidences, in fusion 
as well as in deep- inelastic reactions [H, IS @> 13 we have 
now a good understanding of the process and stimulating 
new perspectives from the use of radioactive beams. 

During the charge equilibration process taking place in 
the first stages of dissipative reactions between colliding 
ions with different N/Z ratios, a large amplitude dipole 
collective motion develops in the composite dinuclear sys- 
tem, the so-called dynamical dipole mode. This collec- 
tive dipole gives rise to a prompt 7-ray emission which 
depends: i) on the absolute value of the intial dipole mo- 
ment 
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being Rz and Rn the center of mass of protons and of 
neutrons respectively, while R p and Rt are the projectile 
and target radii; ii) on the fusion/deep-inelastic dynam- 
ics; iii) on the symmetry term, below saturation, that is 
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FIG. 1: The time evolution of the dipole mode in r— space 
D{t) (solid lines) and p~ space DK(t) (dashed lines, in /m _1 ) 
and the correlation DK(t) — D(t) at incident energy of 
6AMeV and 9AMeV for b = 2/ra. 



acting as a restoring force. 

A detailed description is obtained in a microscopic 
approach based on semiclassical transport equations, of 
Vlasov type, introduced in the nuclear dynamics in col- 
laboration with David |10| . where mean field and two- 
body collisions are treated in a self-consistent way, see 
details in ^lj. Realistic effective interactions of Skyrme 
type are used. The numerical accuracy of the transport 
code has been largely improved in order to have reliable 
results also at low e nerg ies, just above the threshold for 
fusion reactions 0, Il3| The resulting physical picture 
is in good agreement with quantum Time-Dependcnt- 
Hartree-Fock calculation [lij . In particular we can study 
in detail how a collective dipole oscillation develops in 
the entrance channel [T^| . 

First, during the approaching phase, the two partners, 
overcoming the Coulomb barrier, still keep their own re- 
sponse. Then follows a dinuclear phase where the relative 
motion energy, due to the nucleon exchange, is converted 
in thermal motion and in the collective energy of the 
dinuclear mean field. In fact the composite system is not 
fully equilibrated and manifests, as a whole, a large am- 
plitude dipole collective motion. Finally thermally equi- 
librated reaction products are formed, with consequent 
statistical particle/radiation emissions. 

We present here some results for the 32 S + 100 Mo 
(N/Z asymmetric) reaction at 6 and 9 AMeV, recently 
studied vs. the "symmetric" 36 S + 96 Mo counterpart in 
ref.0. In FigH (left columns) we draw the time evolu- 
tion for b — 2/m of the dipole moment in the r-space 
(solid lines), D(t) = ^£-X(t) and in p— space (dashed 

lines), DK(t) = 11/ ft, where II = - with 

Pp (Pn) center of mass in momentum space for protons 
(neutrons), is just the canonically conjugate momentum 
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FIG. 2: Density plots of the neck dynamics for the 32 S + 100 Mo 
system at incident energy of GAMeV and 9AMeV. 



of the X coordinate, see 0, 0|. On the right 
hand side columns we show the corresponding correla- 
tion DK(t) — D{t) in the phase space. We choose the 
origin of time at the beginning of the dinuclear phase. 
The nice "spiral-correlation" clearly denotes the collec- 
tive nature of the mode. From FigQ we note that the 
" spiral-correlation" starts when the initial dipole moment 
D(t = 0), the geometrical value at the touching point, is 
already largely quenched. This is the reason why the 
dinucleus dipole yield is not simply given by the " static" 
estimation but the reaction dynamics has a large influ- 
ence on it. 

A clear energy dependence of the dynamical dipole 
mode is evidenced with a net increase when we pass from 
6AMeV to 9AMeV. A possible explanation of this ef- 
fect is due to the fact that at lower energy, just above 
the Coulomb barrier, a longer transition to a dinuclear 
configuration is required which hinders the isovector col- 
lective response. From Fig|21 a slower dynamics of the 
neck during the first 40/m/c — 60/m/c from the touch- 
ing configuration is observed at 6 AMeV. When the 
collective dipole response sets in the charge is already 
partially equilibrated via random nucleon exchange. 

The bremsstrahlung spectra shown in FigEI support 
this interpretation. 

In fact from the dipole evolution given from the 
Vlasov transport we can directly apply a bremsstrahlung 
("bremss" ) approach ^jjj to estimate the "direct" photon 
emission probability (E-y — Hco): 
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where D"(ui) is the Fourier transform of the dipole accel- 
eration D"(t). We remark that in this way it is possible 
to evaluate, in absolute values, the corresponding pre- 
equilibrium photon emission. In the same FigOwe show 
statistical GDR emissions from the final excited residue. 
We see that at the higher energy the prompt emission 
represents a large fraction of the total dipole radiation. 

In the Table we report the present status of the Dy- 
namical Dipole data, obtained from fusion reactions.. We 



3 



TABLE I: The percent increase of the intensity in the linearized -y-ray spectra at the compound nucleus GDR energy region 
(the energy integration limits are given in the parenthesis), the compound nucleus excitation energy, the initial dipole moment 
D(t = 0) and the initial mass asymmetry A for each reaction. 



Reaction 


Increase (%) 


E*(MeV) 


D{t = 0) (fm) 


A 


Ref 


4U Ca+ iUU Mo 


16 (8,18) 


71 


22.1 


0.15 




36 S+ I04p d 




71 


0.5 


0.17 




lb O+ MS Mo 


36 (8,20) 


110 


8.4 


0.29 
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48 Ti +6 4 Ni 




110 


5.2 


0.05 




^S+ itJU Mo 


1.6 ±2.0 (8,21) 


117 


18.2 


0.19 




36 S+ 96 Mo 




117 


1.7 


0.16 


^S+ luu Mo 


25 (8,21) 
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FIG. 3: The bremsstrahlung spectra for the 32 S + 100 Mo sys- 
tem at incident energy of QAMeV and 9AMeV (solid line) 
and the first step statistical spectrum (dashed line) for three 
impact parameters. 



note the dependence of the extra strength on the in- 
terplay between initial dipole moment and initial mass 
asymmetry: this clearly indicates the relevance of the 
fusion dynamics. 

We must add a couple of comments of interest for the 
experimental selection of the Dynamical Dipole: i) The 
centroid is always shifted to lower energies (large defor- 
mation of the dinucleus); ii) A clear angular anisotropy 
should be present since the prompt mode has a definite 
axis of oscillation (on the reaction plane) at variance with 
the statistical GDR. At higher beam energies we expect 
a further decrease of the direct dipole radiation for two 
main reasons both due to the increasing importance of 
hard NN collisions: i) a larger fast nucleon emission that 
will equilibrate the isospin before the collective dipole 
starts up; ii) a larger damping of the collective mode. 
This has been observed in ref.|8| and more exps. are 
planned 

Before closing I would like to note two interesting 



developments for future experiments with radioactive 
beams: 

• The prompt dipole radiation represents a nice cool- 
ing mechanism on the fusion path. It could be a 
way to pass from a warm to a cold fusion in the 
synthesis of heavy elements with a noticeable in- 
crease of the survival probability, [T^ . 

• The use of unstable neutron rich projectiles would 
largely increase the effect allowing a detailed study 
of the symmetry potential, below saturation, re- 
sponsible of the restoring force of the dipole oscil- 
lation El 



III. ISOSPIN PHYSICS IN A COVARIANT 
APPROACH 

We move now to a relativistic framework for the de- 
scription of the isovector part of the effective interac- 
tion. I will focus then the attention on relativistic heavy 
ion collisions, that provide a unique terrestrial opportu- 
nity to probe the in-medium nuclear interaction in high 
density and high momentum regions. An effective La- 
grangian approach to the hadron interacting system is ex- 
tended to the isospin degree of freedom: within the same 
frame equilibrium properties (EoS, [l9|) and transport 
dynamics (pol l2l|) can be consistently derived. Within 
a covariant picture of the nuclear mean field, for the de- 
scription of the symmetry energy at saturation (04 pa- 
rameter of the Weizsaecker mass formula) (in a sense 
equivalent to the a± parameter for the iso-scalar part), 
extracted in the range from 28 to 36 MeV, there are dif- 
ferent possibilities: (a) considering only the Lorentz vec- 
tor p mesonic field, and (b) both, the vector p (repulsive) 
and scalar S (attractive) effective fields [22J, |23j, |24| . The 
latter corresponds to the two strong effective cu (repul- 
sive) and a (attractive) mesonic fields of the iso-scalar 
sector. We get a transparent form |2ll23|: 
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with E* F = y fc| + m* 2 . 

Once the 04 empirical value is fixed from the p — 5 
balance, important effects at supra-normal densities ap- 
pear due to the introduction of the effective S field. In 
fact the presence of an isovector scalar field is increasing 
the repulsive p-meson contribution at high baryon densi- 
ties 0, via a pure relativistic mechanism, due to the 
different Lorentz properties of these fields (the vector p 
field grows with baryon density whereas the scalar 5 field 
is suppressed by the scalar density). Dynamical, non- 
equilibrium, effects can be more sensitive to such "fine 
structure" of the isovector interaction. We will see the 
vector couplings give 7-boosted Lorentz forces, and so 
we expect a larger isospin dependence of the high energy 
nucleon propagation. Moreover, the scalar S field natu- 
rally leads to an effective (Di rac) mass splitting between 
protons and neutrons YA \TA 121 IM l2aj2fil. with influ- 
ence on nucleon emissions and flows [23, UM- In order 
to explore the symmetry energy at supra-normal densi- 
ties one has to select signals directly emitted from the 
early non-equilibrium high density stage of the heavy ion 
collision. A transverse momentum analysis is important 
in order to select the high density source [2^, |3(| ■ The 
description of the mean field is important, since nucleons 
and resonances are dressed by the self-energies. This will 
directly affect the energy balance (threshold and phase 
space) of the inelastic channels. 



(upper signs for neutrons), where p s — p sp + p sn , j a = 

3p + jnr Ps3 = Psp - Psn, j" = 3p ~ 3n are the total and 

isospin scalar densities and currents and f a ,u,p,s are the 
coupling constants of the various mesonic fields. cr(p s ) 
is the solution of the non linear equation for the a field 

□E3 

For the description of heavy ion collisions we solve 
the covariant transport equation of the Boltzmann type 
|2fi EH ] within the Relativistic Landau Vlasov(i?Ly) 
method, phase-space Gaussian test particles |3lj. and 
applying a Monte-Carlo procedure for the hard hadron 
collisions. The collision term includes elastic and inelas- 
tic processes involving the production/absorption of the 
A(1232MeF) and iV*(1440MeV}resonances as well as 
their decays into pion channels, |32ll33[. 

A relativistic kinetic equation can be obtained from 
nucleon Wigner Function dynamics derived from the ef- 
fective Dirac equation [2l|. The neutron/proton Wigner 
functions are expanded in terms of components with def- 
inite transformation properties. Consistently with the 
effective fields included in our minimal model one can 
limit the expansion to scalar and vector parts F^ (x, p) = 

Fg\x,p) + j fl F^ fl (x,p), i = n,p. We get after some 
algebra a transport equation of Vlasov type for the 
scalar part ffap**) = F^/M*; {p^d* + \p* vl F^ v + 

M*(<9^M*)]cV*}/;(2;,P* M ) = 0, with the field tensors 
Ff v = d»pf + d v p* >1 . The trajectories of test particles 
obey to the following equation of motion: 



Relativistic Transport 

The starting point is a simple phenomenological ver- 
sion of the Non-Linear (with respect to the iso-scalar, 
Lorentz scalar a field) Walecka model which corresponds 
to the Hartree or Relativistic Mean Field (RMF ) ap prox- 
imation within the Quantum-Hadro-Dynamics [l9(. Ac- 
cording to this model the baryons (protons and neutrons) 
are described by an effective Dirac equation ("f^k*^ — 
M*)^(x) = 0, whereas the mesons, which generate the 
classical mean field, are characterized by correspond- 
ing covariant equations of motion. The presence of the 
hadronic medium modifies the masses and momenta of 
the hadrons, i.e. M* = M + S s (effective masses), 
k*^ = k^ — S M (kinetic momenta), where we have in- 
troduced the scalar and vector self-energies E s , E M . 
For asymmetric matter the self-energies are different for 
protons and neutrons, depending on the isovector me- 
son contributions. We will call the corresponding models 
as NLp and NLpS, respectively, and just NL the case 
without isovector interactions. We will show also some 
results with Density Dependent couplings, in order to 
probe effects which go beyond the RMF picture. For 
the more general NLpS case the self-energies of protons 
and neutrons read: 



E a (p,n) = -f a <r{ps) ± fsPs3, 
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d 

Tr 1 



M*(x) 
M*(x) 



Fr(x t (r))+&*M;(x) . (5) 



In order to have an idea of the dynamical effects of 
the covariant nature of the interacting fields, we write 
down, with some approximations, the "force" acting on 
a particle. Since we are interested in isospin contributions 
we will take into account only the isovector part of the 
interaction [34j |: 



dpi 
dr 



Piv 



Jp M* 



VJ% - d v J 3 



±/p^Vp 3 T/ 5 Vp53, (p/n) (6) 



The Lorentz force (first term of Eq.JBJ) shows a 7 = 
boosting of the vector coupling, while from the second 
term we expect a 7-quenched S contribution. We remark 
that the Lorentz- like force is absent in the non-relativistic 
Vlasov transport equation discussed before. This nicely 
shows the qualitative different dynamics of a fully rela- 
tivistic approach. You cannot get it just inserting a rel- 
ativistic kinematics in the classical transport equations. 
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evident is the splitting in the high p t region of the elliptic 
flow. In the (p + S) dynamics the high-p t neutrons show 
a much larger squeeze — out. This is fully consistent with 
an early emission (more spectator shadowing) due to the 
larger repulsive p- field. The V 2 observable, which is a 
good chronometer of the reaction dynamics, appears to 
be particularly sensitive to the Lorentz structure of the 
effective interaction. We expect similar effects, even en- 
hanced, from the measurements of differential flows for 
i light isobars, like 3 H vs. 3 He. 



FIG. 4: Differential neutron-proton flows for the 132 Sn + 124 
Sn reaction at 1.5 AGeV (b=6fm) from the three different 
models for the isovector mean fields. Top: Transverse Flows. 
Bottom: Elliptic Flows. Full circles and solid line: NLpS. 
Open circles and dashed line: NLp. 



IV. COLLECTIVE FLOWS 

The flow observables can be seen respectively as the 
first and second coefficients of a Fourier expansion of the 
azimuthal distribution [3q: 4¥(y,Pt) ~ 1 + 2Vi cos{(j>) + 



2V2Cos(2(j)) where p t — yp^+p y is the transverse mo- 
mentum and y the rapidity along beam direction. The 
transverse flow can be also expressed as: Vi(y,p t ) = {—)■ 
The sideward (transverse) flow is a deflection of forwards 
and backwards moving particles, within the reaction 
plane. The second coefficient of the expansion defines 

-). It measures 



"if 



the elliptic flow given by V2(y,pt) 1 
the competition between in-plane and out-of-plane emis- 
sions. The sign of V2 indicates the azimuthal anisotropy 
of emission: particles can be preferentially emitted ei- 
ther in the reaction plan e ( V2 > 0) or out-of-plane 
(squeeze — out, V% < 0) [35l l3t| . The pt-dependence 
of V2 is very sensitive to the high density behavior of the 
EoS since highly energetic particles (p t > 0.5) originate 
from the initial compressed and out-of-equilibrium phase 
of the collision. For the isospin effects the neutron-proton 

differential flows V^~ p '(y,p t ) = V{ l 2 (y,pt)~Vf 2 (y,p t ) 
have been suggested as very useful probes of the isovec- 
tor part of the EoS since they appear rather insensitive 
to the isoscalar potential and to the in medium nuclear 
cross sections, |37ll3^| . 

In heavy-ion collisions around lAGeV with radioac- 
tive beams, due to the large counterstreaming nuclear 
currents, one may exploit the different Lorentz nature of 
a scalar and a vector field, see the different 7-boosting 
in the local force, Eq.JSJ). In Fig 0] transverse and el- 
liptic differential flows are shown for the 132 Sn + 124 Sn 
reaction at 1.5 AGeV (b = 6/to), that likely could be 
studied with the new planned radioactive beam facilities 
at intermediate energies, [34|. The effect of the different 
structure of the isovector channel is clear. Particularly 



V. ISOSPIN EFFECTS ON SUB-THRESHOLD 
KAON PRODUCTION AT INTERMEDIATE 
ENERGIES 

Particle production represent a useful tool to constrain 
the poorly known high density behaviour of the nuclear 
equation of state (EoS) [33, ■ In particular pion and 
(subthreshold) kaon productions have been extensively 
investigated leading to the conclusion of a soft behaviour 
of the EoS at high densities, Kaons (K°,K+) 

appear as particularly sensitive probes since they are pro- 
duced in the high density phase almost without subse- 
quent reabsorption effects |4j. At variance, antikaons 
(K°,K~) are strongly coupled to the hadronic medium 
through strangeness exchange reactions 0, . Here we 
show that the isospin dependence of the K° :+ production 
can be also used to probe the isovector part of the EoS: 
we propose the K°/K + yield ratio as a good observable 
to constrain the high density behavior of the symmetry 
energy, E sym , 

Using our RMF transport model we analyze pion and 
kaon production in central 197 Au + 197 Au collisions in 
the 0.8 — 1.8 AGeV beam energy range, with different 
effective field choices for E sym . We will compare results 
of three Lagrangians with constant nucleon-meson cou- 
plings (NL... type, see before) and one with density de- 
pendent couplings (DDF, see [24|), recently suggested 
for better nucleonic properties of neutron stars |45| . In 
the DDF model the f p is exponentially decreasing with 
density resulting in a rather "soft" symmetry term at 
high density. In order to isolate the sensitivity to the 
isovector components we use models showing the same 
"soft" EoS for symmetric matter. 

Pions are produced via the decay of the A(1232) reso- 
nance and can contribute to the kaon yield through colli- 
sions with baryons: irB — > YK. All these processes are 
treated within a relativistic transport model including an 
hadron mean field propagation. The latter point, which 
goes beyond the "collision cascade" picture, is essential 
for particle production yields since it directly affects the 
energy balance of the inelastic channels. 

Fig. [3 reports the temporal evolution of A ±,0,++ res- 
onances and pions (n ,0 ) and kaons (K + '°) for central 
Au+Au collisions at 1 AGeV. It is clear that, while the 
pion yield freezes out at times of the order of 50/m/c, 
i.e. at the final stage of the reaction (and at low densi- 
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FIG. 5: Time evolution of the A ±,0,++ resonances and pi- 
ons tt '° (left), and kaons K°' + (right) for a central (6 = 
fm impact parameter). Au+Au collision at 1 AGeV incident 
energy. Transport calculation using the NL, NLp, NLpS and 
DDF models for the iso-vector part of the nuclear EoS are 
shown. 



ties) , kaon production occur within the very early stage of 
the reaction, and the yield saturates at around 20 fm/c. 
Kaons are then suitable to probe the high density phase 
of nuclear matter. From Fig. [SJwe see that the pion re- 
sults are weakly dependent on the isospin part of the nu- 
clear mean field. However, a slight increase (decrease) in 
the 7T~ (tt + ) multiplicity is observed when going from the 
NL (or DDF) to the NLp and then to the NLpS model, 
i.e. increasing the vector contribution f p in the isovector 
channel. This trend is more pronounced for kaons, see 
the right panel, due to the high density selection of the 
source and the proximity to the production threshold. 
The results for the DDF model, density dependent cou- 
plings with a large f p decrease at high density, are fully 
consistent. They are always closer to the NL case (with- 
out isovector interactions) but the difference, still seen 
for 7r + ' _ , is completely disappearing for K ' + , selectively 
produced at high densities. 

When isovector fields are included the symmetry po- 
tential energy in neutron-rich matter is repulsive for neu- 
trons and attractive for protons. In a HIC this leads 
to a fast, pre-equilibrium, emission of neutrons. Such a 
mean field mechanism, often referred to as isospin frac- 
tionation Q, Q , is responsible for a reduction of the neu- 
tron to proton ratio during the high density phase, with 
direct consequences on particle production in inelastic 
NN collisions. 

Threshold effects represent a more subtle point. The 
energy conservation in a hadron collision in general has 
to be formulated in terms of the canonical momenta, i.e. 
for a reaction 1 + 2 — » 3 + 4 as Sj„ = (k^ + k% ) 2 = 
(&3 + &4 ) 2 = s out . Since hadrons are propagating with 
effective (kinetic) momenta and masses, an equivalent 
relation should be formulated starting from the effective 
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FIG. 6: 7r~ /tt + (upper) and K + /K° (lower) ratios as a func- 
tion of the incident energy for the same reaction and models 
as in Fig. In addition we present, for Ei, earn — 1 AGeV, 
NLp results with a density dependent p-coupling (triangles), 
see text. The open symbols at 1.2 AGeV show the corre- 
sponding results for a 132 Sn + 124 Sn collision, more neutron 
rich. 



in- medium quantities k*^ = £+ — S M and m* = m + E s , 
where £ s and £ M are the scalar and vector self-energies, 
Eqs.Q). In reactions where nucleon resonances, espe- 
cially the different isospin states of the A resonance, and 
hyperons enter, also their self energies are relevant for 
energy conservation. We specify them in the usual way 
according to the light quark content and with appro- 
priate Clebsch-Gordon coefficients ^jj^. The self-energy 
contribution to the energy conservation in inelastic chan- 
nels will influence the particle production in two differ- 
ent ways. On one hand it will directly determine the 
thresholds and thus the multiplicities of a certain type of 
particles, in particular of the sub-threshold ones, as here 
for the kaons. Secondly it will modify the phase space 
available in the final channel. 

In fact in neutron-rich systems mean field and thresh- 
old effects are acting in opposite directions on particle 
production and might compensate each other. As an ex- 
ample, nn collisions excite A _, ° resonances which decay 
mainly to 7r~. In a neutron-rich matter the mean field 
effect pushes out neutrons making the matter more sym- 
metric and thus decreasing the n~ yield. The threshold 
effect on the other hand is increasing the rate of 7r~'s 
due to the enhanced production of the A~ resonances: 
now the nn — > pA~ process is favored (with respect to 
pp — > nA ++ ) since more effectively a neutron is con- 
verted into a proton. Such interplay between the two 
mechanisms cannot be fully included in a non-relativistic 
dynamics, in particular in calculations where the baryon 
symmetry potential is treated classically ^Tj ■ 

In the 0.8 — 1.8 AGeV range the sensitivity is larger 
for the K°/K + compared to the ir~ /tt + ratio, as we can 



see from FigEJ This is due to the subthreshold produc- 
tion and to the fact that isospin effects enter twice in 
the two-steps production of kaons, . Between the two 
extreme DDF and NLpS isovector interaction models, 
the variations in the ratios are of the order of 14 — 16% 
for kaons, while they reduce to about 8 — 10% for pions. 
Interestingly the Iso-EoS effect for pions is increasing at 
lower energies, when approaching the production thresh- 
old. We have to note that in a previous study of kaon 
production in excited nuclear matter the dependence of 
the K°/K + yield ratio on the effective isovector interac- 
tion appears much larger (see Fig. 8 of ref.J!^). The point 
is that in the non-equilibrium case of a heavy ion collision 
the asymmetry of the source where kaons are produced 
is in fact reduced by the n — > p "transformation" , due 
to the favored nn — > p/Sr processes. This effect is al- 
most absent at equilibrium due to the inverse transitions, 
see Fig. 3 of ref. 33]. Moreover in infinite nuclear matter 
even the fast neutron emission is not present. This result 
clearly shows that chemical equilibrium models can lead 
to uncorrect results when used for transient states of an 
open system. In the same Fig. [5] we also report results 
at 1.2 AGeV for the 132 Sn + 124 Sn reaction, induced 
by a radioactive beam, with an overall larger asymmetry 
(open symbols). The isospin effects are clearly enhanced. 



VI. TESTING DECONFINEMENT AT HIGH 
ISOSPIN DENSITY 
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The hadronic matter is expected to undergo a phase 
transition into a deconfined phase of quarks and gluons 
at large densities and/or high temperatures. On very 
general grounds, the transition's critical densities are ex- 
pected to depend on the isospin of the system, but no ex- 
perimental tests of this dependence have been performed 
so far. Moreover, up to now, data on the phase transi- 
tion have been extracted from ultrarelativistic collisions, 
when large temperatures but low baryon densities are 
reached. In order to check the possibility of observing 
some precursor signals of some new physics even in col- 
lisions of stable nuclei at intermediate energies we have 
performed some event simulations for the collision of very 
heavy, neutron-rich, elements. We have chosen the reac- 
tion 238 U + 238 U (average proton fraction Z/A = 0.39) 
at 1 AGeV and semicentral impact parameter b — 7 fm 
just to increase the neutron excess in the interacting re- 
gion. To evaluate the degree of local equilibration and 
the corresponding temperature we have followed the mo- 
mentum distribution in a space cell located in the cm. 
of the system; in the same cell we report the maximum 
mass density evolution. Results are shown in Fig. [7| We 
see that after about 10 fm/c a nice local equilibration is 
achieved. We have a unique Fermi distribution and from 
a simple fit we can evaluate the local temperature. At 
this beam energy the maximum density coincides with 
the thermalization, then the system is quickly cooling 
while expanding. In Fig0 lower panel, we report the 



FIG. 7: Uranium-Uranium 1 AGeV semicentral. Correlation 
between density, temperature, momentum thermalization in- 
side a cubic cell 2.5 fm wide, located in the center of mass of 
the system. Lower panel: density, temperature, energy den- 
sity, momentum, proton fraction. Curves in the upper-left: 
black dots - baryon density in po units; grey dots - quadrupole 
moment in momentum space; squares - resonance density. 



time evolution of all physics parameters inside the cm. 
cell in the interaction region.. We note that a rather ex- 
otic nuclear matter is formed in a transient time of the or- 
der of 10 fm/c, with baryon density around 3 — 4po, tem- 
perature 50 — 60 MeV, energy density 500 MeV fm~ 3 
and proton fraction between 0.35 and 0.40, well inside 
the estimated mixed phase region, see the following.. 

A study of the isospin dependence of the transition 
densities has been performed up to now, to our knowl- 
edge, only by Mueller |4^. The conclusion is that, mov- 
ing from symmetric nuclei to nuclei having Z/A ~ 0.3, 
the critical density is reduced by roughly 10%. Here we 
explore in a more systematic way the model parameters 
and we estimate the possibility of forming a mixed-phase 
of quarks and hadrons in experiments at energies of the 
order of a few GeV per nucleon. 

Concerning the hadronic phase, we use the rela- 
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tivistic non-linear Walecka-type model of Glendenning- 
Moszkowski (GM...) where the isovector part is 

treated just with p-meson couplings, and the iso-stiffer 
NLpS interaction |5l). 

For the quark phase we consider the MIT bag model 
|52j with various bag pressure constants. In particular we 
are interested in those parameter sets which would allow 
the existence of quark stars |53U54j. i.e. parameters sets 
for which the so-called Witten-Bodmer hypothesis is sat- 
isfied [55ll56| . One of the aim of our work it to show that 
if quark stars are indeed possible, it is then very likely 
to find signals of the formation of a mixed quark-hadron 
phase in intermediate-energy heavy- ion experiments |5l| . 

The scenario we would like to explore corresponds to 
the situation realized in experiments at moderate energy, 
in which the temperature of the system is at maximum of 
the order of a few ten MeV. In this situation, only a tiny 
amount of strangeness can be produced and therefore we 
can only study the deconfincmcnt transition from hadron 
matter into up and down quark matter. Since there no 
time for weak decays the environment is rather different 
from the neutron star case. 




0.1 0.2 0.3 0.4 o.e 



FIG. 8: Variation of the transition density with proton frac- 
tion for various hadronic EoS parameterizations. Dotted line: 
GM3 parametrization; dashed line: NLp parametrization; 
solid line: NLpS parametrization. For the quark EoS, the 
MIT bag model with B 1/4 =150 MeV. The points represent 
the path followed in the interaction zone during a semi-central 
i32 gn+ i32 Sn collision at i AGeV (circles) and at 300 AMeV 
(crosses), see text. 



Mixed phase 

The structure of the mixed phase is obtained by impos- 
ing the Gibbs conditions \5Q, |6(| for chemical potentials 
and pressure and by requiring the conservation of the 
total baryon and isospin densities 

(H) _ (Q) (H) _ (Q) 

Pb — Pb ) A*3 — Ma ) 
^)(T,^)=F(«(T,,g), 

PB = (1 - X)Pb + XPb i 

P3 = (1-X)P?+XP?, (7) 

where \ is the fraction of quark matter in the mixed 
phase. In this way we get the binodal surface which 
gives the phase coexistence region in the (T, Pb, Pz) space 
For a fixed value of the conserved charge p 3 , 
related to the proton fraction Z/A = (1 + p 3 /pb)/2, we 
will study the boundaries of the mixed phase region in the 
(T, ps) plane. We are particularly interested in the lower 
baryon density border, i.e. the critical/transition den- 
sity p cr , in order to check the possibility of reaching such 
(T, p cr , p 3 ) conditions in a transient state during an HIC 
at relativistic energies. In the hadronic phase the charge 
chemical potential is given by /13 = 2E sym (pB)^ ■ Thus, 
we expect critical densities rather sensitive to the isovec- 
tor channel in the hadronic EoS. 

We compare the predictions on the transition to a de- 
confined phase of the two effective Lagrangians GM 3 and 
NLpS. The isoscalar part is very similar while the isovec- 
tor EoS is different, because in GM3 we only have the 
coupling to the vector p-field. In Fig. [S] we show the 
crossing density p cr separating nuclear matter from the 
quark-nucleon mixed phase, as a function of the proton 



fraction Z/A. We can see the effect of the ^-coupling 
towards an earlier crossing due to the larger symmetry 
repulsion at high baryon densities. The ^-exchange po- 
tential provides an extra isospin repulsion of the hadron 
EoS, and its effect shows up in a further reduction of the 
critical density. 

In the same figure we report the paths in the (p, Z/A) 
plane followed in the cm. region during the collision of 
the n-rich 132 Sn+ 132 Sn system, at different energies. We 
see that already at 300 AMeV we are reaching the border 
of the mixed phase, and we are well inside it at 1 AGeV. 



Deconfinement Precursors 

Statistical fluctuations could help in reducing the den- 
sity at which drops of quark matter form. The reason is 
that a small bubble can be energetically favored if it con- 
tains quarks whose Z /A ratio is smaller than the average 
value of the surrounding region. This is again due to the 
strong Z/A dependence of the free energy, which favors 
clusters having a small electric charge. These configura- 
tions can easily transform into a bubble of quarks having 
the same flavor content of the original hadrons, even if 
the density of the system is not large enough to allow 
deconfinement in the absence of statistical fluctuations. 
Moreover, since fluctuations favor the formation of bub- 
bles having a smaller Z/A, neutron emission from the cen- 
tral collision area should be suppressed, what could give 
origin to specific signatures of the mechanism described 
in this paper. This corresponds to a neutron trapping ef- 
fect, supported also by a symmetry energy difference in 
the two phases. In fact while in the hadron phase we have 
a large neutron potential repulsion (in particular in the 
NLpS case), in the quark phase we only have the much 
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smaller kinetic contribution. If in a pure hadronic phase 
neutrons are quickly emitted or "transformed" in pro- 
tons by inelastic collisions, when the mixed phase starts 
forming, neutrons are kept in the interacting system up to 
the subsequent hadronization in the expansion stage [5l| . 
Observables related to such neutron "trapping" could be 
an inversion in the trend of the formation of neutron rich 
fragments and/or of the ir~ /ir + , K°/K + yield ratios for 
reaction products coming from high density regions, i.e. 
with large transverse momenta. In general we would ex- 
pect a modification of the rapidity distribution of the 
emitted "isospin" , with an enhancement at mid-rapidity 
joint to large event by event fluctuations. A more de- 
tailed analysis is clearly needed. 



mixed hadron-quark matter at high baryon density in the 
collision, central or semi-central, of neutron-rich heavy 
ions in the energy range of a few GeV per nucleon. As 
signatures we suggest to look at observables particularly 
sensitive to the expected different isospin content of the 
two phases, which leads to a neutron trapping in the 
quark clusters. The isospin structure of hadrons pro- 
duced at high transverse momentum should be a good 
indicator of the effect. 

Many new ideas for fundamental experiments with ra- 
dioactive beams are emerging. My picture of David Brink 
like a tree always starting new braches is more and more 
confirmed 



VII. PERSPECTIVES 

We have shown in few examples the richness of the 
physics we can describe using mean field transport equa- 
tions, inspired by the pioneering works of David M. 
Brink. We have seen that collisions of n-rich heavy ions 
from low to intermediate energies can bring new infor- 
mation on the isovector part of the in-medium interac- 
tion in different regions of high baryon densities, quali- 
tatively different from equilibrium EoS properties. We 
have presented quantitative results for charge equilibra- 
tion in fusion/deep inelastic reactions, differential col- 
lective flows and yields of charged pion and kaon ratios. 
Important non-equilibrium effects for particle production 
are stressed. Finally our study supports the possibility 
of observing precursor signals of the phase transition to a 
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